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Abstract 

In this paper wc examine the IR inflation scenario using the DBI action, where we have 
N muhiple branes located near the tip of a warped geometry. At large N the solutions 
are similar in form to the more traditional single brane models, however we find that it is 
difficult to simultaneously satisfy the WMAP bounds on the scalar amplitude and the scalar 
spectral index. We go on to examine two new solutions where N ~ 2 and N = 3 respectively, 
which both have highly non- linear actions. The sound speed in both cases is dramatically 
different from previous works, and for the = 3 case it can actually be zero. We show that 
inffation is possible in both frameworks, and find that the scalar spectral index is bounded 
from above by unity. The level of non-gaussian fiuctuations are smaller in the N = 2 case 
compared to the single brane models, whilst those in the A^ = 3 case are much larger. 
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1 Introduction 



Although inflation is generally accepted as a model for the early universe, it is still a paradigm 
lacking an explicit derivation from a more fundamental theory such as string theory. Conversely 
we expect that inflationary cosmology provides an important testing ground for the inflationary 
models from string theory thanks to the abundance of precision data [1]. In recent years we 
have witnessed a massive increase in the number of string-cosmology papers (see [2] and the 
references therein). 

One of the interesting and purely stringy models is the Dirac-Born-Infeld (DBI) inflation [3] 
in type IIB string theory. This relies on the existence of a mobile D3-brane moving along the 
throat geometry of a warped flux compactification. The inflaton arises from an open string mode 
representing the relative distance of the mobile brane with respect to some fixed background. The 
non-linear form of the DBI action in this warped background was shown to lead to substantial 
inflation despite the fact that the brane could be moving at near relativistic velocity. It was 
further shown that there are generically two models of this kind of inflation which were dubbed 
the UV and IR models [4]H In the UV model, the brane moves towards an anti-brane located 
near the bottom of a warped throat. The anti-brane would be localised near the bottom of the 
throat since it screens the RR flux present in the compactification, and this corresponds to a 
form of large field inflation. The IR model, in contrast, arises after brane flux annihilation where 
mobile branes travelling down the throat annihilate quantum mechanically with the trapped 
flux creating new branes in the process [6]. These branes would then feel the attractive force of 
branes/fluxes in other throats and so would generically be attracted towards the gluing region. 
This corresponds to a small field inflation model. In fact the most general case is an interpolating 
one, where one considers the brane starting in the UV, travelling down the throat, annihilating 
with the flux and then residual brane returning back up the throat. It was shown that the level 
of non-gaussianities predicted in this model are proportional to the square of the relativistic 
velocity factor. Since this is usually greater than unity in these kind of models, we can falsify 
models of this type by flnding only small non-gaussianities in future experiments. 

It seems unlikely that only a single residual brane emerges from the flux annihilation process, 
since this would require extreme fine tuning of the background fiuxes. More generally one would 
assume that there will be branes left after the flux annihilation, and they all feel the same 
attractive force driving them back up the throat. Indeed this was considered in [4]. However 
here it is supposed that the branes are all separated by distances greater than the string length, 
which just gives rise to a modified version of assisted infiation [7] . Since the branes are all created 
in the same annihilation process, it is natural that they are located at various distances, with 
some being coincident. For simplicity we consider that all the branes are basically coincident at 
the bottom of the throat. We can then describe their dynamics by the non-Abelian DBI action, 
which is given by the Myers action [8] This action is complicated but simplifies when one takes 
the large N limit. Although this appears to introduce large dimensionless parameters into the 
theory, it can still be useful for infiationary model building [6,10]. 

^See also [5] for an interesting alternative. 

''There is also a proposal by Tseytlin [9] for the non-Abelian action, but we use the Myers formalism in this 
paper. 
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Another motivation for considering brane configuration is in relation to reheating [fl]. In 
the usual single brane models, it is assumed that the Standard Model (SM) lives in a separate 
throat from the inflationary throat. It is supposed that once inflation ends, reheating occurs in 
the SM throat through graviton production which tunnels through the internal space. However 
in the absence of brane annihilation at the end of inflation, there is still a residual massless U{1) 
degree of freedom in the inflationary throat. Then the reheating energy can be transferred to 
the U{1) gauge boson rather than to the SM. By contrast, with multiple branes all the reheating 
energy couple directly to the SM degrees of freedom. This is an interesting possibility which 
deserves further investigation, but will be left for future study. 

Since the DBI model of inflation is entirely stringy in origin, it is in our interest to scan as 
much of the possible parameter space of solutions as possible. This will tell us more about the 
robustness and predictiveness of the model, and also whether this kind of inflation is generic. 
There has been recent work in this direction, however we hope to probe more of the parameter 
space by considering multiple branes as they are typically more common in string theory. 

In section 2 we introduce the non-Abelian action and the relevant backgrounds. In section 
3 we consider the simplest approximation of the the large N limit for the action. This shares 
many of the features of single brane models while retaining its inherently non-Abelian structure. 
We investigate the inflationary parameter space and determine the level of the curvature and 
gravitational perturbations. In section 3, we discuss the case of finite N. To highlight the 
differences between multi and single brane inflation, we examine the cases of A = 2 and A = 3 
in some detail. We calculate the speed of sound in both cases, and use this to constrain the 
parameter space. We also calculate the relevant inflationary parameters and the perturbation 
amplitudes and non-gaussianities of the model. We finish with some conclusions and suggestions 
for future work. 

The results of this paper can be summarized as follows. 

In the context of the non-Abelian world volume theory of A D3-branes, we find, in the large 
A limit, similar solutions to those of the single brane scenarios previously considered, with the 
additional constraint that the background charges must be larger than the number of branes. 
Our model differs significantly in some respects from a single probe brane, despite sharing many 
similarities. The main difference is that now the world-volume theory is playing the role of the 
universe, rather than just the inflaton sector in the standard IR approach to infiation. In fact 
this makes the model a hybrid between that of DBI infiation [3] and Mirage Cosmology [12]. 

With certain assumptions, we have found that infiation is reasonably generic, i.e. we can 
generate at least 60 e-folds with relative ease due to large A suppression. We investigate in 
detail the cases of Ads and mass-gap backgrounds and find that the level of non-gaussianities 
are the same and that spectral indices are positive running and bounded from above by unity. 
This large A^ model may, however, be ruled out by experiment as that the large factors of A" 
tend to enhance the size of scalar perturbations. 

The case of finite A is particularly interesting and striking. For A = 2 we find that the 
non-gaussianities are suppressed a little relative to the single brane case. The extra stringy 
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degrees of freedom present in the theory serve to induce a so-called 'fuzzy' potential term which 
makes the form of the action for finite N significantly different from that of a single brane. As a 
consequence we find that in the relativistic approximation, the speed of sound becomes a rather 
complicated function at finite A'^ and we are restricted to certain regions of the moduli space of 
solutions. For the N = 3 case we even find that the speed of sound actually has zeros in its 
function. The physical explanation for this is not immediately obvious to us, but is clearly a 
signature of the underlying complexity of the non-Abelian world volume theory. 



2 The non-Abelian action 



In this section, we introduce the non-Abelian action [8], which provides an effective description 
of coincident D-branes. We first discuss the string theory background relevant for our model. 

Let us assume that the type IIB spacetime metric factorises as follows 

ds"^ = h^dsl + h-'^{dp^ + ds\^), (2.1) 

where the four-dimensional metric is taken to be the usual Friedmann-Robertson- Walker (FRW) 
form characterised by the scale factor o(t), a throat region is over some five-dimensional manifold 
X5, and h \s a, warp factor. We take the so-called IR inflation, where the branes are initially 
localised at the IR tip of the throat where the warp factor is small. These IR localised branes 
arise naturally in the context of brane-flux annihilation. 

Now consider a background where we have M units of D3-brane flux threading some three- 
cycle of the internal space. If we insert N' l)3-branes into this background such that they fill the 
large 3+1 dimensional spacetime, they will feel an attractive force from the flux and roll down the 
throat to annihilate quantum mechanically. Then there will be precisely N = M — N' D3-branes 
created after this annihilation process. This provides the initial conditions for our model. In 
realistic string compactifications, there are many throats, and their respective orientifold images. 
Therefore it is natural that these remnant branes will feel attractive force from anti-branes in 
other throats. Thus the remnant branes will travel along the throat towards the internal space. 

This model has been investigated recently for a single brane, N = 1. However this means that 
we need fine tuning to ensure that there is only a single brane remaining after flux annihilation. 
In a more generic situation, several branes will remain after this process. For simplicity, we 
consider the case in which the branes are at distances smaller than the string scale with U{N) 
symmetry. Then we should use the non-Abelian DBI effective action which describes coincident 
branes. 

The bosonic components of the Myers action [8]: 

S = -T3 j d^eSTr {j-Aet{Eab + Ea^{Q~^ - Sy^Ejb + XFab) ^detQ]^ . (2.2) 
The fields Eat are the non-Abelian pullback of the linear combination of closed string fields 
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Eab = Gab + Bab, while the matrix Q is determined as 

Qi = Si+iX[cl)\<P']Ekj, (2.3) 

with the 0* being scalar fields on the world-volume of the D-branes corresponding to the trans- 
verse fluctuations and A = 27r/^ with Ig the fundamental string length. Implicitly the scalars 
are now N x N matrices transforming in the irreducible representation of their respective gauge 
group. This irreducibility condition essentially means that they are in their lowest energy state. 
The symmetrised trace prescription means that we must trace over the symmetric average of all 
the fields in the action. Building upon recent work [16, 21] we assume that the bulk B fields is 
zero near the tip of the throat, and consider the case where the transverse coordinates define a 
fuzzy (S*^ embedded in a three cycle in the X5 manifold. 

A fuzzy sphere is defined in a similar way to an ordinary sphere, except that the defining 
relation for the d-dimensional sphere 

d+l 

^(x^)2 = R^, x' E ]R'^+\ (2.4) 

i=l 

is subject to the replacement of x* by A^-dimensional matrices x*. Any function expanded in 
terms of spherical harmonics on the sphere is subject to a cut-off due to the finite size of the 
matrices. Mathematically this means that we are truncating the algebra of functions on the 
sphere. Thus a fuzzy sphere is essentially a sphere where we can identify only N points. For the 
two-sphere, it is easy to see that in the large limit, the fuzzy S*^ coincides with the classical 
two-sphere. Note that this is not generally the case for higher dimensional fuzzy spheres. 

Now that we have restricted ourselves to the S0{3) ^ SU{2) algebra, we expect the scalars 
to be proportional to its generators. As argued in [8], it is preferable for the generators to be 
in the irreducible representation as this corresponds to the lowest energy configuration. The 
matrix ansatz that we take for our scalar fields is thus = Rd^ [i = 1,2,3), where it! is a 
variable of canonical mass dimension, and the a* are the A^-dimensional generators of the SU{2) 
algebra. 

In flat space, we can find solutions where the radius of the fuzzy sphere grows without bound. 
However for most compactifications there will be a maximum bound on the size of the fuzzy 
sphere, as illustrated by the famous example of the AdSs x metric where the fuzzy sphere 
radius is bounded by the radius of the five-sphere. In the warped backgrounds where the length 
of the throat provides the cutoff scale, the throat is smoothly glued onto the five-dimensional 
internal space. We call this the UV cut-off, analogous to the UV brane in the Randall-Sundrum 
models [13]. 

Let us insert the metric (12. ID into the Myers action ()2.2p . After expanding all the determi- 
nants, we find that the effective action for coincident D3-branes in this backgrounc^ becomes 

5 = -Ta y d^eSTr (h'^aWl - h-^X'^a^a^R'^Vl + ^X'^a^a^h-^R'^ - a^/i^lN + a^ViR)!^^ , 

(2.5) 

^This is based upon earlier work [14]. 
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where the second term arises as the leading order contribution from the Chern-Simons couphng 
of the bulk RR fields, while the final term is a flux induced potential which comes from fluxes 
present in the background. Both terms are singlets under the trace, which is why they appear 
multiplied by the N x N identity matrix. The potential generated will depend on the topology 
of the internal space, and also the quantisation constraints of the fluxes. Here we have absorbed 
a factor of Tg"^ into the potential to make it dimensionless, and set the dilaton to be unity. This 
agrees with the supergravity background generated by D3-branes, and also the tip solution of 
the Klebanov-Strassler throat [15,16]. 



We see that the non-linear form of the DBI action leads to an infinite series expansion in 
powers of the SU{2) generators. The action is thus complicated by the symmetrised trace over 
the generators. The action can be simplified by keeping only the leading term in each of the 
square roots, which is the large limit. In this approximation, all the terms in the action are 
proportional to the identity matrix and Tr(a*a*) = CTr(l) where C is the quadratic Casimir of 
the A^-dimensional representation of SU{2). 



We now minimally couple the DBI action to the standard Einstein-Hilbert action. Fortu- 
nately the Einstein part of the action arises naturally in the process of compactification, so we 
can consistently include it in our model. It is convenient to introduce the notation 

W{R, C) = Vl + 4X^h-^a^a'R^, (2.6) 

which is essentially an additional potential induced by the fuzzy sphere geometry, which we 
will refer to as the fuzzy potential. For a single probe brane, this contribution is always unity, 
since the corresponding matrix representation is commutative. We can now determine the non- 
zero components of the energy-momentum tensor on the world-volume. It gives the energy and 
pressure densities as 

P = -TgSTr (h'^W{R, C)Vl - h-^X^a^a^R? - /i^ + . (2.8) 

These form the basis for all our analysis in the subsequent sections and we will return to them 
in due course. 



One of the reasons that we can have DBI inflation in this background is that the speed of 
sound is usually very small compared to unity. In the standard canonically normalised slow 
roll models, this factor is always unity since the scalar field moves slowly. The speed for the 
non-Abelian case is 

C2 dP/dk 
" dE/dk' ■ 
which is valid provided that the entropy corrections are negligible. 
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3 The large N limit 



In this section, we focus on the large N limit. We show that the solution shares many of the 
features with the single brane models. 

Let us restrict ourselves to the large N limit when we ignore the backreactior^. The large N 
limit has proven to be useful for inflation in many other contexts [6,10]. Now we can approximate 
the symmetrised trace by a trace once we neglect the contributions of terms of the order 
in the action. The resulting action simplifies to 

S = -T3 j (fiNa^ [h^^ 1 - h-^X^CR^^ 1 + 4A2C'/i-4i?4 _ /j4 + V{R)^ (3.1) 

where is the tension of the DS-branes given by 

with Ms being the mass scale for the open strings and gs being the asymptotic string coupling 
which we take to be small to allow for perturbatively defined strings. In our paper we will assume 
that the coupling is set to Qs ~ 10~^ in order to make order of magnitude approximations. The 
back reaction effect will be small provided that we ensure MK >> N is also satisfied, in addition 
to the large assumption. It is instructive to make a redefinition of the scalar field. Firstly 
we switch to 'physical' coordinates using the relation = r^/(A^C') which parameterises the 
physical radius of the fuzzy sphere. This defines the square of the physical radius. Let us also 
define a scalar field (p = r^/T^ with canonical mass dimensions in order to make contact with 
the rest of the literature. 

We have the standard relationship between the four-dimensional Planck scale and the ten 
dimensional one through the volume of the warped space Vq: 

Ml = \. (3.3) 

Interestingly it was shown in [18] that with minimal assumptions about the volume of the throat, 
one can find the following bound on the maximal field variation: 

where MK is the contribution from the background fluxes. In our case we demand that MK >> 
N to neglect the back reaction upon the geometry, which restricts our field to move over very 
small distances in Planckian units. 



This amounts to a constraint on the energy of the D-branes which could lead to an unphysical solution upon 
compactification. However it can be interpreted as a metric constraint on the bulk fluxes. 
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We find from (12.20 the equation of motion of the </> field: 

= w{(t>,c)^^4> + ?,H^w{4),c)^ + ^^'^^ 



Ti\'^CW{(t),C) V h 
+ 2nh'W{^, (7)7 (^2h^ - - mh^h' + nV' , (3.5) 

where primes are derivatives with respect to (j) and we have also introduced 

7= ; ^ =, (3.6) 

1 - h~^T.^^(j)^ 

for the analogue of the relativistic factor for the DBI action. This implies that the velocity of 
the brane is bounded as 

02 < h^Ts. (3.7) 

These last two equations are exactly the same as in the case of a single brane [3]. Recall that 
when we take the large limit of the fuzzy S'^, we recover the classical S'^ with units of 
charge. 

Using the general expression (j2.9p . we calculate the speed of sound = 1/7^, in agreement 
with that of single brane inflation. As in the single brane model, we now assume that the scalar 
field is monotonic, at least for the early times. This assumption allows us to switch to the 
Hamilton-Jacobi formalism, whereby we use the field itself as the dynamical parameter rather 
than time variable. We differentiate the Friedmann equation = E/3Mp with respect to 
time, dropping terms proportional to (j), and use the continuity equation, E = —3H{P + E) to 
get 

^ - " (3.8) 



N-iW{(l)) ' 



where the fuzzy potential W is now an explicit function of (j), and H' is the derivative of the 
Hubble parameter with respect to the inflaton. Substituting this into p.6p . we obtain 



We can use (j3.9p to write the velocity of the inflaton without reference to the relativistic factor 
7: 

-2M^H' 

^ (3.10) 



Ar2^2 + 4M4i7'2/i-4T3-i 



which allows us to write the corresponding Priedmann equation solely as a function of the 
inflaton: 

= ^{W{cP)h\cPh{c^)+V{^)-h\cP)) (3.11) 
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which is reminiscent of the equation found in [3] for the Abehan DBI inflation model. The main 
difference here is the presence of factors of N and the fuzzy sphere induced potential W{(j)), the 
latter being an inherently non-Abelian feature. In this manner we have written the variables 
purely as functions of the scalar field. 



We may be concerned that the DBI models of inflation do not exhibit standard attractor 
solutions for inflation, since we expect relativistic motion. To check this, let us suppose that 
i?o(0) is a solution of ()3.1ip . which can be either inflationary or non-inflationary. We add to 
this a linear homogeneous perturbation 6H{cj)). The attractor condition will be satisfied if it 
becomes small as (j) increases. Upon substituting H = Hq + SH into (13. lip and linearizing the 
resultant expression, we find that the perturbation obeys 



H',5H' = ^J^H,5H, (3.12) 

which has the general solution 



2M2 



5H{(j)) = 5H{(f>i) exp 



(3.13) 



where 6H(4>i) is the value at some initial point (pi, and 7 = j{Hq). Because Hq and d(j) have 
opposite signs, the integrand within the exponential term is negative definite, and all linear 
perturbations indeed die away. This means that there is an attractor solution for this model 
regardless of the initial conditions and the velocity of the brane. This is also true for the single 
brane solutions [4] with = 1 = W{(j)). 



We note that the equation of state in this model is drastically different from the canonical 
field models. It is for the large N 



If the potential dominates all the other terms, we recover the usual de-Sitter solution with 
w ~ — 1. However the DBI admits more interesting solutions due to its non-linear nature. For 
example, if we consider ultra-relativistic motion where 7 >> 1, and demand that the /i^ terms 
are suppressed, we obtain 

which can be very small depending on the scale of the fuzzy potential, and may give rise to 
a matter phase in the asymptotic velocity limit. This shows that we have a larger parameter 
space of solutions for lv than in the standard inflationary scenarios. 



Our analysis thus far has been general. To make more detailed investigation of the infla- 
tionary signature of this model, we must determine the background potential. Let us consider 
nonzero fluxes inducing the warped throat solution. The coincident branes localised at the bot- 
tom of the throat will feel a potential generated by branes/fluxes in other throats and will move 
towards them. We then expect a tachyonic potential of the form 




(3.16) 
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with higher order even powers of (p because of the Z2 symmetry of the throat. The various 
constants wiU be determined by the choice of compactification, and also non-perturbative effects. 
The IR DBI inflation is thus a special case of small field inflation. The constant Vq is the scale set 
by the fluxes, and need be large to be able to neglect back-reactive effects in our model. This is 
significantly different from the UV inflationary model. In this paper we are mainly interested in 
the IR solution, and we consider that most of the dynamics will take place in a region dominated 
by the potential energy. 



3.1 Inflationary observables and constraints 

In this subsection, we focus on inflationary solutions in two specific backgrounds. We show that 
inflation can be achieved in this model by allowing an appropriate tuning of the parameters in 
the theory. 

We wish to consider inflation near the tip of a warped throat where (j) is small. However 
we are also interested in solutions where the background is no longer approximately AdS but 
is constant. This is motivated by the work of [18] and also the fact that a finite warp factor 
appears to be a generic feature of the IR end of warped throats. This is a phenomenologically 
motivated solution, because once we compactify our theory on a compact space, we expect the 
fiuxes to back-react on the bulk geometry forming a throat. We expect that these throats will 
be of the Klebanov-Strassler (KS) variety, which have a finite cut-off at the origin. This cut-off 
is generically exponentially small due to its dependence on the three-form fluxes. 

In order to mimic a constant warping in our non-compact theory we choose to put in 'by 
hand' a constant warp factor parameterised by a mass term ^, where we expect that fx is small 
at a scale set by the bulk fluxes: 

M^) = (3.17) 

where we have used L to denote the background charge. It should be noted that this L will be 
different from that in the AdS-like backgrounds. When (p goes to zero, the warp factor remains 
finite. Strictly speaking away from the origin there may be a different (p dependence, but we use 
this form of the warp factor in the following sections. 

The solutions we consider are 

(i) the AdS type cases (where we set /i to 0), or 

(ii) the mass gap cases (where we assume ~ 0). 

Let us assume that infiation occurs very close to the tip of the warped throat, in which case 
we expect the energy density to be dominated by the constant piece of the background potential. 
The other terms are suppressed by the square of the warp factor and will be small in this limit. 
It may appear that the warp factor can be easily vanishingly small, but more care is required 
in cases where h reduces to a constant, since in those backgrounds (as we shall see) the other 
parameters can both be large. Assuming Vq >> h'^{W{(/))'y{(/)) — 1) is satisfied, the Friedmann 
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equation ()3.1ip can be approximated as 



2 NT^Vj^) 

We then find that our solutions for the inflaton velocity and the gamma factor reduce to 



M„V' / T3 



V 3NV ' 
I M2 / y'2 \ 

where we have expressed everything in terms of derivatives of the potential. Clearly these 
expressions are background dependent as they depend on the warp factor h. 



3.2 Inflation in AdS type backgrounds 



Let us consider the locally AdS type backgrounds, such as those studied in the original IR 
inflation scenario [4]. The harmonic function can be approximated hy h ~ (p/L in the near 
horizon region where L is essentially the charge of the background geometry given bj0 = 
gsM KT^ X'^ir'^ /Yolp^^), with M and K the corresponding quanta of flux and Vol(X5) correspond- 
ing to the dimensionless volume of the compact space. We expect to find that Vol(X5) = avr'^, 
where a is a topological parameter. For example we know that a = 1 for the five-sphere and 
a = 16/27 for the manifold T^'^ [18]. In this situation the fuzzy sphere induced potential be- 
comes constant which greatly simplifies the analysis. In fact the second term inside the square 
root is proportional to the ratio of the background fluxes and the number of coincident branes, 
and thus it is not clear a priori whether this term will be small or large. Under our assumption 
of no back reaction, we require the flux term to dominate over N, and so we may expect that 
W >> I which translates into the flux condition 

L2 jv 



Ml » 8^5^ ^ ^/M}fS»iV (3.20) 

in the large N limit. For small values of the string coupling constant, we see that this requires 
the fluxes to be very large. This is to expected from out heuristic arguments regarding the 
backreaction. In the converse limit where W ^ Iwe see that the constraint becomes \^g^MK « 
N - which although provides a bound on is much harder to satisfy within the remit of our 
approximation . 

This requirement also affects the definition of 7: 



^Note we have rescaled this quantity to ensure it has the correct dimensions. This corresponds to the radius 
of curvature for the AdS space scaled by the square of the brane tension 
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As we are interested in DBI inflation, we should take jitj)) >> 1 for the speed of sound to 
be substantially reduced, and so the right hand side will dominate the previous expression. 
Furthermore since we are assuming that the background potential should dominate the energy 
density at this stage of the evolution, we should also assume that Vq >> V2, which implies that 
the new constant piece of 7 will be subdominant. 



In this regime, we can approximately solve the inflaton equation of motion. In fact there is 
a cancellation between terms in the equation which implies that (p ~ 0{(j)^) + 0{(f)^) and so for 
consistency we must drop all terms higher than quadratic in the fields. The resultant expression 
for the field is actually of the same functional form as the single brane case. Actually this is 
expected since we can view the fuzzy sphere as a classical sphere with N units of flux in the 
large limit: 

^ ^ ^, (3.22) 

where the field is initially located at (j) = (f)Q a,t t = t^. This expression on the denominator 
will generally be much smaller than unity as we are assuming that /iq << 1, which is the initial 
value of the warp factor. 



Let us now compute the inflationary parameters in this large N limit. The non-linear form 
of the DBI action prevents us from using the traditional slow roll variables, and so we must 
establish new 'fast roll' variables. In reality the name 'fast roll' is somewhat of a misnomer 
because despite moving relativistically, the non-linear nature of the action allows for the brane 
to be held up on the potential for a significant amount of time as in slow roll scenarios. This 
has already been extensively discussed for the single brane models but the non-Abelian action 
requires us to modify these expressions. Suppose that the leading order term for the epsilon 
parameter expansion is given by 

- = H'^{l-e), (3.23) 



which yields the usual slow roll constraint e = —H/H"^. However now that we are working in 
the Hamilton-Jacobi formalism, we need derivatives with respect to the infiaton. This leads to 
the following modified expressions for the relevant slow roll parameter: 



2M2 /h' 



- 



N-iW{(l3,C) \ H 



(3.24) 



This is clearly equivalent to the usual single brane slow roll conditions where we would have 
N = l,W{(f),C) = 1, although we cannot take the A — > 1 limit in our non-Abelian DBI 
description. Note that this slow roll parameter is suppressed not only by a factor of I/7 as 
in the single brane infiation, but also by an additional factor of 1/A. Intuitively we may have 
expected this since the coincident branes will tend to accelerate much more slowly than a single 
brane. Hence we would expect inflation to last for a longer period of time. 

Using our approximate solutions we flnd for the AdS type backgrounds (and assuming 7 >> 

1) 

^3Mp^ 
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where we have neglected cf/' terms. We must ensure that e < 1 for inflation to occur at all. We 
can also calculate the number of e-foldings for this model using 



~ X^TT / #72 1 - ^ + (^ir) ' (3-27) 



iVe = y Hdt. (3.26) 

We obtain 

where the integration should be between cpQ and the latter being determined though the fast 
roll parameter. Clearly for small (p the first term will dominate the integral and so we drop the 
higher order terms. The result is that the field value A'^e e-folds before the end of inflation can 
be written as 

.... ^/i^^ (3.28) 



L^^/NVo + VS(|)fNeMp' 

which we can use to determine the perturbation spectrum. For completeness, we write the fast 
roll parameter as a function of the number of e-foldings: 



3.3 Inflation in mass gap backgrounds 



The equation of motion for the inflaton in the AdS background is basically the same as in the 
single brane case. Is this also true for the mass gap solution? In this instance the fuzzy sphere 
potential now has non-trivial field dependence which complicates the analysis. Without loss of 
generality we will take the warp factor of the form h ~ fi/L: 



7mg(</') 



1 + ^- 3.30 



In this paper we are mainly interested in the relativistic limits of the theory, and so we 
require 7 to be large. This immediately imposes a constraint on the fuzzy potential - which is 
now an explicit function of the inflaton (j) 

The analytic form of W tells us that it is bounded from below by unity, but has no upper 
bound. Of course in reality we expect W{(j)) to be a monotonically increasing function, however 
for analytical simplicity we will consider the two limits separately. 
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In the first instance let us assume that W ^ 1. Prom (I3.3ip this imphes that we have an 
upper bound on the number of branes 

which can be satisfied by having a small enough warp factor, and suitable hierarchy between 
terms in the potential. Note we are implicity assuming that (j) is non-zero here - since a vanishing 
field corresponds to a fuzzy sphere of zero size due to the relationship between the canonically 
normalised scalar and the fuzzy sphere radius. However we must also impose the constraint 
coming from the definition of the fuzzy potential which actually implies a lower bound on N 
through the relation 

A' » Jgr (3.33) 
SO upon combining these constraints we see that both assumptions are valid provided that 



We now ask about the constraints arising from the converse limit, when both 7 and W are 
large. For the fuzzy potential to be large we must ensure that 

however the relativistic limit also requires to be bounded from below 

7 X 10- Wo 



Combining both of these constraints also results in ()3.34p . It is straight- forward to note that 
due to the constant warp factor near the tip of the throat, the relativistic limit implies that (j) is 
a linear function of time. This was already established for the single brane case by Kecskemeti 
et al [17], and follows from the fact that our assumption requires us to drop the factor of unity 
in the definition of 7, which then uniquely fixes the form of <j) to be constant. The overall scale is 
simply set by the warp factor and the brane tension, and is independent or our parameterisation 
of the fuzzy potential. 

Let us now consider how inflation occurs in this model. Our starting point will once again 
be (I3.24p . Inserting the mass gap solution and also demanding relativistic motion is enough to 
ensure that W drops out of the analysis, thus the solution is independent of the fuzzy potential 
contribution. After integrating to find the number of e-folds we see that the slow roll parameter 
can be written in the form (Nf. e-folds before the end of inflation) 

3NeV2M^h'^ 

We have seen that for a certain range of parameters, we can obtain the required number of 
e-foldings. The real signature of the model lies in the perturbation spectrum and the scale of 
the spectral indices, which we now address. 
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3.4 Cosmological Signatures. 



The derivation of the various perturbation spectra for this model is presented in the Appendix, 
and we refer the reader there to see how the expressions arise. We will simply quote the important 
results in the following section. The main equations we need to calculate the perturbation spectra 
are (|6.13|) and (|6.15|) respectively. In terms of our standard notation employed in the rest of the 
paper these translate into the following conditions when we substitute for the velocity equation 
in the Hamilton-Jacobi formalism, and re-insert the factors of the reduced Planck mass 



A 



2 



-^T ~ 8 (^7^^ ) . (3.38) 



27rMp 



The slow-rolling inflation generally predicts very low non-gaussianity since in the leading order 
the quantum fluctuation are generated by free fields in the dS background. However in the DBI 
inflation, much larger non-gaussianity can be generated since the causality constraint in the 
kinetic term introduces non-linear interactions among different momentum modes of the scalar 
field. Recently, it was shown that in the equilateral triangle limit, the leading-order contribution 
to the non-linearity parameter is given by [18] 



where 



It can be shown that the second term on the right-hand side of (|3.39p vanishes in the large N 
case, as it does in the single brane case. Therefore we have the usual non-gaussian parameter 

0.3272. (3.41) 

This once again emphasizes the similarity between the = 1 and the A^ >> 1 descriptions. 
Current measurements indicate a rather weak bound on the level of non-gaussianities fxL ^ 100, 
however the upcoming Planck mission aims to increase the sensitivity to probe down to regions 
where /tvl < 5. 

We will now discuss the expected level of perturbations in each of the two cases we have 
discussed so far. 



• Ads backgrounds. Using the formulae derived in the appendix we find that the scalar 
amplitude at leading order becomes 
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Now using the WMAP normalisation for the non-gaussianities we must ensure that 7 < 
10\/3, which in turn fixes e through the (weak) constraint that r < 0.5. Therefore we see 
that at horizon crossing we must ensure that 

5^3 WHl^ 

e < -— MKNV2 < (3.43) 

- 16 ^ ^ 

is satisfied where we assume that the string coupHng is roughly 10~^. This essentially 
means that the slow roll parameter must be less than one half. Using this constraint to 
fix the scalar amplitude we find the following condition needs to be satisfied 

Clearly this is linear in NVq, therefore requires that the string scale must be low in order 
for the WMAP normalisation to hold, assuming that the constant part of the potential is 
sub-Planckian. In addition we see that the scalar index can be expanded as a power series 
in e, which at leading order becomes 

~ 1 - 4e - Ae^/^dus + ... (3.45) 

where we have defined the 'perturbation' 



6ns ~ ^ = . (3.46) 



( 102 Mp y/^ 

Recall that the WMAP bound for this index is given by n, = 0.987l:J]:[!;^^. Using the 
constraint in (13.430 we find that the bn^ term must be bounded from below by unity. 
Moreover we see that by constraining e to be smaller - this in fact makes the bn^ term 
larger. Given this, one sees that the scalar index will always be large and negative in this 
instance and therefore incompatible with observation. 

We also see using the relation 

X = S,r (3.47) 



din A; dN, 



e 



evaluated at horizon crossing, that both tensor and scalar indices are positive running. 

Mass Gap Backgrounds. The form of the Hubble parameter in this instance is the 
same as in the Ads case, therefore we expect similar arguments to hold in this instance. 
Once again this favours a smaller string scale using (|3.44|) . Assuming similar constraints 
on r and 7 we see that the only real distinction between the two cases arises through the 
spectral indices, since we again have to ensure that (j3.43p is satisfied. Explicit calculation 
of the scalar index in this case reveals that 

If we restrict ourselves to the case where Vl^ ~ 1 then this simplifies down to the following 
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The scalar index will clearly be sensitive to the magnitude of the second term, which we 
write as —dus/e, and we see that 



> J_ ( 1 _ ^ ) (3.50) 

using the WMAP normalisation. Clearly for small values of the warp factor we will see 
that the 5ns term will be large, which implies that we cannot obtain the observed spectrum 
of scalar perturbations. The only way we can satisfy the experimental data in this instance 
is to assume a large warp factor - however this is contrary to all our assumptions thus far. 
So we must conclude that this particular model does not agree with the data. 

Conversely in the limit where we take W » 1 we find a cancellation between dangerous 
terms which gives the final result for the scalar tilt 

~ 1 - 7£, (3.51) 

which can be seen to arise as the limiting case of the ~ 1 solution. Using the WMAP 
bound, this serves to fix a much tighter bound on e at horizon crossing, which we can 
interpret as a upper bound on 

T/zlO^M^/i^ 

N < (3.52) 

and therefore we can satisfy the experimental bounds on inflation. In both cases we see 
that the indices are positive definite, as in the Ads models. 



3.5 Examples of other solutions 



Let us consider more general solutions which arise from the master equation (|3.11|) . We will 
consider both UV and IR inflationary solutions for generality and then make some comments 
about the general signals of inflation in the large N limit. 



To distinguish between each solution branch, we note that on the UV side we have large 
field inflation. So there will generally be positive contributions to the potential, and the vacuum 
energy Vq will be set to zero. In the IR branch, we have Vq 7^ at all times and the potential 
will generally be taken to be tachyonic indicating that the branes move away from the tip of the 
throat. In the fuzzy sphere picture, the UV side corresponds to a collapsing sphere, while the 
IR side corresponds to an expanding sphere. 



An interesting solution was analysed in [3] when there is no quadratic term in the potential 
ie V{cl)) = Vo — Vi(j)'^. Solving the master equation (|3.11|) for H{cj)) and integrating back to 
find the time dependence of the inflaton, we obtain the following solutions for the AdS type 
backgrounds: 



j,4 



m ^=\i^. (3.53) 
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where the terms in the Hubble parameter are calculated to be 



H^^Jm. H.^m^(}^.V.l (3.54, 



3M2 ' ^ V 3Vb 2Mp V 

Note that there is a potential sign ambiguity in the definition of H4. This is because we can 
consider either the IR inflation (where we have a minus sign in the potential) or the UV inflation 
(where we have a plus sign). In both cases, the equation of motion for the inflaton is the same. 
In the IR case this is the early time, small 4> solution, whereas in the UV case it corresponds to 
the late time solution where the fuzzy sphere has collapsed to almost zero size. In both cases, 
the form of the equation of motion implies that in the small fleld limit the inflaton is moving 
relativistically with large 7. The number of e-foldings can be determined as follows: 

iVe ~ Hoite - to) + {4>\te) " ^\to)) , (3.55) 

where to and te are the times at the beginning and end of inflation, respectively. We must ensure 
that tf > te, where tf represents the time at which we can no longer trust our approximations. 

Let us now reconsider the mass gap solution. Since the throat is flnite and the warp factor 
is actually constant for small values of (p, there is no reason why the field is moving rapidly near 
the origin. Let us take a quadratic potential of the form V{(l)) = Vq i V2(j)'^, where the positive 
sign signifies UV type infiation. For small values of the inflaton, we have the Hubble parameter 
i?((/)) = Hq + H2(p'^, where the coefficients are determined in a similar way to those in the AdS 
case. If we assume that the field is small (in both cases), the general solution can be written as 



This solution is valid for the IR inflation. For the UV solution, one only need substitute a minus 
sign in front of the V2 term inside the square root. Again this means that for the IR solution 
this is the early time evolution, while for the UV solution it is the late time evolution of the 
Hubble parameter. In both cases, the solution for the field can be written as 

4M^H2(tf-t) 

<^(t)~(/>oe N , tf>t, (3.57) 
which means that the field is rolling non-relativistically because H2 and (j) are both small. 



Let us focus initially on the IR solution. If we wish the inflaton to be increasing, correspond- 
ing to the branes moving away from the origin, we are forced to choose the minus sign in H2. 
The other sign is a solution where the fleld is getting smaller. In the event that V2 is zero, we 
flnd either a solution where the inflaton is at a constant value, or H2 is positive deflnite and 
again the fleld is getting smaller. The former case corresponds to de-Sitter expansion since we 
can immediately integrate the solution for the scale factor to get a = e^"*. In the UV region 
we must also impose an additional reality constraint 9Hq > 8NT3V2. If this reality bound is 
saturated, the field is again rolling towards the origin. In both cases the evolution of the field 
is essentially determined by the vacuum energy Vq which sets the overall scale for the Hubble 
parameter. 
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Using the solution to the equation of motion, we can calculate the number of e-foldings and 
write is as a function of the field. In the UV case the branes would be near the bottom of the 
throat, and so this correspond to late time evolution of the inflaton. As such we interpret this 
as the early stage evolution of the IR inflationary model. We obtain 

iVe ~ Ho{te - io) + ^ (<A'(ie) " ^'(io)) , (3.58) 

e-foldings in this regime and we must again ensure that tg ^ tf- Since we expect (j) to be 
small over this time period in accordance with our approximation of the Hubble parameter, the 
dominant contribution to the number of e-foldings comes from the constant part of the potential. 
However since N » 1 we may well see a sizable contribution to the number of e-foldings. 



Let us consider what happens in more generality, although we will assume that the Hubble 
parameter is generically of the form shown in ()3.18p . Focusing our attention on the power 
spectra, we find that the amplitudes are given by 

Al = -oT72> (3-59) 



7r2M2' 



_ g2 / m^H'^ 



1/2 



3 



16e 1 + 



4Mj^g^2 x-i/2 



where we have included r as the ratio of the tensor /curvature amplitudes. Recall that each of 
these is to be evaluated at horizon crossing if we wish to normalise them to the WMAP data [1]. 
Note that e is expected to be small at horizon crossing. We can repeat the same analysis as 
before and consider limits of the term in parenthesis, (i) In the first case when we consider 
relativistic motion, the scalar amplitude and ratio reduce to 



,2 H^\H' 

At • 



Air'^eNWh'^./T^ 
8eNWh'^y/% 



M^\H'\ 



(3.60) 



If we satisfy the condition that r < 0.24 then we find from the scalar amplitude, saturating the 
bound, that 8H^ ~ 10~^M2 or more concretely that 

iVm~10-^°(^)\ (3.61) 

assuming that the string coupling is around 10^2 _ This will be extremely difficult to satisfy 
under the assumption of large N, and also vacuum energy dominance, unless we are willing to 
postulate a low string scale 
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(ii) If we consider the non-relativistic limit, we find the solutions 



The first expression implies that e » at horizon crossing, which can be satisfied provided 
that the energy density of the inflaton is vanishingly small. In fact if this condition is met, the 
ratio will simultaneously be satisfied. A quick calculation shows that if we demand r < 0.24 (ie 
using the strongest possible WMAP bound), then the energy density must satisfy E < 10~^Mp 
which is very small in Planck units. 

Using the fact that /tvl < 100 we can obtain a bound on N for our theory. For an arbitrary 
warp factor we see that 



To an order of magnitude approximation the numerical factor is 0(1). Now for mass gap 
backgrounds we may generally expect W ^ 1 or W >> 1 and /i is a small constant which forces 
N to be large. This will be further enhanced by a smaller string scale unless the derivative of 
the Hubble parameter is vanishingly small. For the Ads scenarios we have competition between 
the W and /i^ terms in the denominator - so we would expect N to be set explicitly by the 
choice of background fluxes and the string scale. 

The general form for the scalar spectral index can be calculated to give 

ns-1 = -2e(^2-^A^ (3.64) 

for an arbitrary warp factor. In principle the backreaction effects will appear through a redefi- 
nition of the harmonic function, and so this expression should be valid for all theories satisfying 
our assumptions. This equation simplifies once we assume relativistic motion, ie 7 >> 1 

For AdS type solutions the fuzzy potential is constant, however we can see that the limiting 
solutions are 

H" h' 

Clearly in our simplest case analyzed in the previous section we see that for W >> 1 the second 
term will be identically zero thus cancelling out any dangerous dependence. For the mass 
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gap solution we find similar expressions to those presented above except that the h' terms will 
be zero - at least to leading order. Unless the backreaction dramatically alters the solution, we 
should expect inflation to favour the W » 1 regime, in which case the scalar index is essentially 
only a function of the potential and its derivatives - the overall scale being set by the e term. 

In general we expect the number of e-foldings to be enhanced by factors of thus making 
the universe generically very flat. However the factors of tend to increase the level of scalar 
and tensor perturbations, making it difficult to satisfy observational bounds without imposing 
restrictive fine tuning. 



4 Inflation at flnite N 

In this section, we investigate cosmic inflation due to a small number of coincident branes. We 
begin with some general remarks about the finite formalism before specialising to two simple 
cases, namely N = 2 and = 3. In these cases the action is highly non-linear and gives an 
expression for the speed of sound and the inflationary parameters in certain regions of the phase 
space. They are very different from those in the single brane models. 



4.1 General remarks and motivations 

We will now switch to the finite N formulation of the non-Abelian Myers action, using the 
prescription for the symmetrised trace as given in [19]. We believe this prescription to be correct, 
however a concrete proof remains an outstanding problem. In most models of brane inflation, 
the bulk fluxes are tuned so that only a single brane is left after the brane-flux annihilation 
process. In the context of the landscape, this is a very special case and the general expectation 
is that there remain several residual branes which will tend to coincide to minimise their energy 
in a warped throat. In the first part of the paper we looked at the large A^ limit, which has many 
problems due to the large back-reactive effects on the geometry although it is a more general 
solution than the single brane cases. In the remainder of the paper we will look at the solution 
when there are a handful of residual branes at the tip of a warped throat. This means that we 
can effectively neglect back-reactive effects as in the single brane models, while still retaining 
the enhanced non-Abelian world-volume symmetry. 

Let us rewrite the expressions for the energy and pressure of the coincident branes in a more 
suitable manner for finite A^: 
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where we have used the following definitions 



X = \^h-\ Y = 4X^R^h-\ (^)= ^.o/o ^^^{!! V> (4-2) 

\mj r(3/2 -m)r(m + 1)' ^ ^ 

and the fact that the potential is a singlet under the trace. We employ the symmetrisation 
procedure in [19]. The basic formulas we need are then 

{n+l)/2 

STr[(a*a*)'"] = 2(2m + 1) ^ (2i - 1)"* 

i=l 
n/2 

= 2(2m + l)^(2z)'". (4.3) 

i=l 

The first line corresponds to odd n = N — 1, and the second line to even n. Note that we move 
from working with the A^-dimensional representation to the spin representation with n = 2J. 
It is important to consider what we mean by a physical radius in this context. We use the 
definition 

.^ = A^;i^L,n.„_( ^^^' )=AW. (4.4, 

which implies that the Lagrangian will converge for velocities from to 1, and moreover that 
the radius of this convergence will be unity. This definition is consistent with what we know 
about the solution in the large N limit. 



To illustrate the additional complexity arising from the finite N solution, let us calculate the 
speed of sound in two examples using ()2.9p . For = 2, we find 



' 3 + y(4-xi?2) ' ^ ■ ^ 

where = r^/A^. This is obviously far more complicated than the large expression which 
is = 1 - Xii^. For iV = 3, we obtain 

C^fr^r - 3) - - ^Xii')i3 + 16Y - 8XR^[l + 6Y]) 

' 3 + i6y(i - xi?2) ' ^ ■ ^ 

where now 

^2 ^ jsjote that in both cases, we recover the usual result that (7^ = 1 

when the velocity of the branes is zero. 



We may have anticipated the fact that at large N the solution should resemble that of a single 
brane, through knowledge of the Myers effect. However the finite N solutions would be expected 
to be radically different. One could imagine trying to start with the supergravity dual of this 
model, i.e some fivebrane wrapping a two-cycle with N unnits of C/(l) flux. However in order to 
capture finite N corrections to the theory requires the use of the world volume star-product, and 
not simply the ordinary product assumed here. Extrapolating backwards to extremely small 
values of requires us to keep higher order terms in the expansion of the star product, which 
significantly complicates the form of the resulting Lagrangian. 
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4.2 Two brane inflation. 



In this subsection, we consider inflation driven by two coincident branes moving in the warped 
throat. This gives rise to a [/ (2) symmetry on the world-volume. The relevant expressions for 
the energy and pressure can be calculated using (14. ip and (14. 3p : 

E = 2tJ^^1±^^^^^ + V-A (4.7) 

Note that in order to keep the energy finite, we should impose the constraint cp"^ < fi^T^. This 
also ensures that the contribution coming from the DBI part of the action will be non-negative, 
as can be seen from the first term in the numerator of the energy equation. 

In general it is difficult to get solutions due to the complicated form of the energy density. 
So let us make the approximation that the inflaton is rolling ultra relativistically. We can define 
the relativistic factor 7 much as we did in the large N solution by 7 = (1 — XB?)'^/"^ . 

We now write the energy and pressure as functions of 7, and then take the large 7 limit. By 
utilising the conservation equation and dropping all acceleration terms, we can find the solution 

MlH' 

7 ~ T , (4.8) 

where H' is the derivative of the Hubble parameter with respect to (j). The sign in (|4.8p cor- 
responds to the two choices (f) = ±\/T^h? + ... in the expansion of the velocity (f) about its 
saturation value. The — sign is for the choice (/> > while the -|- sign for (/> < 0. Note that in 
order to have 7 > 0, we must demand that H' < for > and H' > for (p < 0. The choice 
of sign here is vital to obtaining the correct solution branch for infiation. 

We can rewrite the speed of sound as a function of 7 and Y as follows: 

7' V 37^(1 + F) + y J ' ^ ^ 

which shows the finite N corrections to the equation in this limit. Without recourse to a specific 
background, we can make the following observation: If we consider limiting solutions for Y, ie 
that it is either >> 1 or << 1, all y dependence drops out of the expression and the equation 
reduces to C| ~ 1/(87^). The sound speed is thus only a third of that in the large N limit. 
This appears to be the attractor point for the velocity with this action regardless of background 
choice. In this case we find (I3.40p is given by 

2{1 -XR'^) {3 + iY -XYR"^) 
In the 7^1 case, the non-linearity parameter is 

/tvl- 0.2472, (4.11) 
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which is a httle smaller than the large (and single brane) solution and effectively means that 
we can satisfy the observational bounds whilst considering larger velocities than in the single 
brane scenario. 



The corresponding Friedmann equation in this case is 

E 

H'^ = ^, (4.12) 

where we are using the energy as defined in ()4.7p . Substitute our expression for 7 into this 
equation, we find that the Hamilton Jacobi equation for = 2 (in the large 7 limit) becomes 

(4.13) 

for an arbitrary flux induced potential. At this stage, we could either specify the form of the 
Hubble parameter and then consider how this modifies the potential [17], or we could specify 
the form of the potential and then solve for H. We use this latter approach as this appears to 
be more within the spirit of the Hamilton-Jacobi formalism we have employed thus far. 



4.2.1 Inflation in AdS type backgrounds 



Let us first consider a solution where the potential is dominated by a constant, V ^ Vq. Let us 
also assume that the background is approximately AdS. In the small field limit, we expect the 
/i^ term is negligible compared to the remaining terms. So as a first approximation, we ignore 
its contribution. Solving this differential equation, with an appropriate constant of integration 
C, we obtain 

^.(,).?||..„,.(g^(_l,C«). (4.14) 

Since H' < 0, we have assumed > in obtaining this solution. In the limit of very small (j)^ we 
can approximate the solution by ~ 2T^Vq / {?>M'^) , which is of the same functional form as 
in the large A^ case. However this expression is not consistent with the our approximation that 
7 >> 1 since the ratio H' /(fp' is approximately zero for vanishingly small values of 4>. Therefore 
we must be careful to choose a regime of validity where this solution is valid. Careful inspection 
shows that the function 7^ has a turning point in the small field limit, with a maximal value 
given by 

^ (4.15) 



(max 



9 L4yo'/'^/^+T' 



where Y is constant in the AdS background, and the value of the field at this point is given by 
(f) = (;/)max ~ L^\/6Vo/(4Mp). Thus to consider inflation in this region, we need some 'extreme 
fine tuning' to set up the initial value of the field. 

Rather than proceeding this way, we make the Taylor series expansion of the Hubble pa- 
rameter for small (j) but without dropping the term in (j4.14p . This means we must include 
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quartic terms in the expansion of the Hubble parameter and so we need 



Hi(l))^^H,<p\ (4.16) 



i=0 



We will also keep quartic terms in the inflationary potential for consistency, V{<j)) ~ Vq — 
V2(l)'^/2 - V^cji^/A + .... Equating the various coefficients, we find that the linear term in the 
Hubble parameter actually vanishes, leaving us with the residual terms 



^0 = (4.17) 



H2 



3M^ 
T3V2 



H4 = — ^ ( T/4i^r.f ^ + 4t!^^^ - UM^H^L^Ts + 6M^L'^^/I^H^ 

l2M^L^Ho,/T3 \ ^ 3 P d J p V d 2 

We find that the constant piece of the potential dominates the Hubble parameter when the 
field is vanishingly small. The sign of the last term is potentially ambiguous which can lead to 
interesting cosmological behaviour. It turns out that the Hubble term is extremised at the usual 
(j) = solution (which is a local maximum), and there exists a non-trivial solution given by 



8H4 



-3Hs ± ^9i7| - 32H2H^ , (4.18) 



where we must require the term inside the square root to be non-negative. We can rewrite this 
reality constraint as 

Ho 8HnL'^ , , 

TT ^ (4-19) 

We now find the possibility of a 'cosmic turnaround' because H2 will be negative definite for 
those regions of phase space where H4 is also negative. If we concentrate on regions of (p near 
the origin, H' < so we necessarily have </) > in order for 7 > in our approximation ()4.8p . 
Of course, we have implicitly assumed that the field is monotonic so we cannot say anything 
about the reality of such a bounce solution within the current framework. 

It is easy to solve the equation of motion in this limit for relativistic motion. As in the case 
of a single brane, and for the large N solution, we obtain the following term for the infiaton 
equation of motion 

, 00 



1 - <PoVn{t - to)/L^ ' 

where again we define 4>o as the field value at time t = to, and it can be seen that (p > 0. 



(4.20) 



To see the implications of this for infiation, we must first determine which are the relevant 
parameters in this finite formulation. The modified 'fast roll' parameter in this case can be 
written as 
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where the sign is related to the sign for (p and we need to demand that inflation ends when e = 1 
as usual. Recall that for small (f), H' < and we need to choose the minus sign in this equation 
(coming from the choice (f> > 0). Inserting our expression into the Hubble parameter, we see 
that £ can be expanded in powers of the inflaton. Keeping only the leading order term (which 
amounts to dropping 0{(j)'^) contributions), we see that inflation will end around 



2 u2 \ 1/3 







2\Ho 



(4.22) 



The corresponding number of e-foldings given by this Hubble parameter is generically a power 
series in (p. We expect the dominant contribution to arise from the constant piece Hq as in the 
standard inflationary scenario. Integrating over the field, we find the expressions for the inflaton 
as a function of e-folding number: 

Inserting this back into the 'fast roll' parameter (I4.2ip . we find that the dependence on the 
number of e-foldings is of the same functional form as in the large N case: 

For the perturbation amplitudes we can use the general results developed in the appendix. We 
find that the gravitational wave amplitude will be constant to leading order, and given by 

the corresponding expression for the scalar amplitude is given by 




We choose to re-write this in terms of the tensor amplitude as 



1/3^ 



3 



^|~10"V?^|l + 60<|^| j (4.27) 

Now the non-gaussianity condition implies that the bound 1 << 7^ << 400 must be satisfied 
to comply with observation. The spectral indices for this model at large 7 are 



~ (4-28) 
4X2 

^ -(1 + 60X,)+--- 

where the first line is understood to be evaluated at horizon crossing, and we have written 
X2 = {2)Mp / L^V2Y^^ a dimensionless parameter. Note that the tensor index is negative, but 
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suppressed by the Hubble parameter. Clearly the scalar index is bounded from above by unity, 
and so normalising to WMAP data implies that < X2 < 0.05, or more concretely that 

^ > M^. (4.29) 



2.4 X 104 



Using this constraint in ()4.27p we see that the term in brackets varies between 1 and 64. The 
relationship between the two amplitudes is characterised by the parameter r and so we recover 
the anticipated DBI relation 

r ~ -. (4.30) 

7 

The fact that 1 << 7 < 20 in this model implies that r will generically be small and thus well 
within the WMAP confidence bounds. Furthermore this implies that the tensor amplitude will 
be smaller in magnitude than the scalar one, something like 10"^'' for a range of 7. Using the 
normalisation for the scalar amplitude, namely that it satisfies ^| ~ 10~^ at horizon crossing, 
we see that this constrains the potential in terms of the string scale. For small X2 we find 



^00(10^) (M,^^ 



gs \Mp 



> 1 (4.31) 



where we are interested in an order of magnitude approximation. Whilst for the maximal value 
of X2 we recover the same constraint but with an additional factor of 10^ in the numerator. 
Clearly both solutions are sensitive to the mass splitting between Planck and string scales and 
imply that generic inflation prefers the string scale to be close to the Planck scale in order 
not to have super-Planckian scalar potentials. For example if we have Ms ~ 10~^Mp then the 
potential constraint becomes 0(10^ — 10^)Vo > 1. What about the constraint in (14.290 ? Upon 
substituting for the background parameters we see that this equation can be written as 

V2 > ^ (4.32) 

- MK M4 ^ ^ 

which sets the inflaton mass scale. 



4.2.2 Inflation in mass gap backgrounds 



If we repeat the analysis for the mass gap backgrounds, (assuming that the constant part of the 
potential dominates), we find the following solution for the Hubble parameter: 



-(*)-/^ v'^^'^j 

where we have used the fact that the warp factor is constant to write the solution as a function 
of /i. 

Substituting this into the gamma constraint, we must require that the solution is larger than 
unity even when (p is vanishingly small. We make a Taylor series expansion of the resultant 
function, and find that 7^ oc sech^(F(C)), where the amplitude of the function is determined 
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by the ratio of the potential and the warp factor. Also note that F{C) is an function of the 
Casimir, whose precise form is not relevant for understanding the physics of the solution. Now 
the hyperbolic trigonometric function is a decreasing function of its argument, which forces us 
to take the limit « h^/Vo in order for the large velocity expansion to hold. 



Let us assume that we can in fact take this limit and consider the implications for inflation. 
Calculation of the fast roll parameter e yields the following 

which is a decreasing function of the inflaton field. Thus after some critical field value (pc we 
will find a solution where infiation never ends. It may appear that this is an artifact due to the 
neglect of higher order terms in the potential. However if we consider quartic terms in V {(/)), 
and also up to the same order in a Taylor expansion of H{(j)), we find the same result that £ 
is a decreasing function of the infiaton. We conclude that in the relativistic limit that inflation 
(once started) never end^ unless we turn on extra effects, such as nonzero gauge fields. This 
result is not anything that we expect from the results of the single brane case [3], or the large A'" 
limit discussed in the previous sections and appears to be a distinctly finite N effect. Of course, 
it may well be that standard infiation can occur for moderate values of 7. However one should 
probably need numerical results to see it. 



Of course, the mass gap background will eventually give way to something similar to the 
Ads solution, where the harmonic function will have explicit dependence upon the infiaton field. 
Therefore we expect inflation to end in this regime. The fact that the mass gap solution has flnite 
warping means that it will be relatively easy to produce the necessary 60 e-folds of expansion. 



4.2.3 Non-relativistic limit 



Let us restrict ourselves to the non-relativistic regime in order to see the consequences for 
brane inflation It would of course be more preferable to obtain an interpolating solution 
between these two extremes, however it is analytically challenging and would be better suited 
to a numerical analysis. After performing a series expansion of the continuity equation, we find 
the following solution for the velocity of the field in the Hamilton-Jacobi formalism: 

^ = - (3 + 4/) ' ^'-^'^ 
which means that the corresponding Friedmann equation reduces to 

Tf; 7^ 1 + 2y + J,^ + V{(j)) - /l^ (4.36) 



^Technically this is no longer true once the branes reach the gluing region. However the effective action is no 
longer expected to be a good description of the physics in this region. 
®This has recently been examined in [20]. 
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where = {1 + Y)/{2> + AY). Let us solve this equation by considering a standard quadratic 
potential. There are two solution branches, one of which has an imaginary component of H. 
We ignore this solution as it appears unphysical. The other real solution can be parameterised 
by a quadratic Hubble parameter with nonzero components given by 



For inflation to occur we must ensure that we take the minus sign in the solution for H2. 
The speed of sound in this instance reduces to 



where we can substitute (p for derivatives of the Hubble parameter. Of course our non-relativistic 
expansion assumes that the subleading term is smaller than unity in order for the speed of sound 
to remain a real function. This constrains the field to satisfy 

(t> « ^^T^ -■ 4.39 

15 

The non-relativistic assumption means that we can find inflationary solutions even in the mass 
gap backgrounds. Using the definition of the e parameter we find that the leading order contri- 
bution yields 

e ~ e^^^^ (4.40) 

at e-folds before the end of inflation. We have introduced the dimensionless ratio j3 = 
MpH2/ Hq for simplicity. For inflation to occur we must ensure < 0, or that H2 < 0. However 
the fact that the slow roll parameter is now exponential implies that the level of scalar per- 
turbations will now be enhanced by this exponential term. Our assumptions assume that (p is 
small, so the speed of sound is essentially unity in this instance. It transpires that the simplest 
equation to study is the tensor to scalar ratio r, which is now given by r ~ 16e. If we demand 
that r < 1/4 to satisfy the more stringent bound, then after some algebra we find the following 
constraint 

p 

where we have explicitly left in the numerical value. If this bound is not satisfied then we find H2 
to be very small which suppresses the number of e-foldings and sufficient infiation is generically 
difficult to achieve. 

For Ads type backgrounds we have the following solution for the Hubble parameter 



where we must again take the minus sign for an infiationary solution. Interestingly for Y « 1 
we see that the solution becomes exactly the same as the mass gap one. We can follow the same 
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procedure and obtain a similar result for the number of e-folds. The difference is of course due 
to the constant nature of the factor Y. 



e ~ exp 



(^) 



where 5 is defined in a similar way to (3 except that the Hubble parameters are different in this 
case. 

The observed constraint on the amplitude ratio can now be written as 

which is slightly different from the mass gap solutions. The function F{Y) acts to increase the 
rhs of the expression above, ranging from F{Y) ~ 3 ^ 4\/y depending upon our choice of 
fluxes. Larger values of Y impose clearly impose tighter constraints on the parameter V2 and so 
one would anticipate that smaller values are more preferential. 

In both instances we see that in order to satisfy the scalar curvature constraints, we require 
the dominant term in the potential to satisfy the following 



which can be combined with the expressions for V2 to yield a constraint purely on that variable 
in terms of the Planck and string scales (and also the fluxes for the Ads case). We have again 
assumed a string coupling of Qs ~ 10~^ in the above expression. 



4.3 Three brane inflation 



Let us now move to the case where there are three coincident branes, giving rise to a U{2>) 
world-volume symmetry. We have the energy and pressure: 

E = 2Ti\ ^ — ^ : — ^ + , (4.46) 

lVrT2F(l-2Xi?2)3/2 2 ' ' ^ ' 



p = -2n 



/i''(i + 4y-4Xi?2[i + 3y]) 

a/i - 2Xi?Vl + 2Y 




The symmetry breaking induced by a gauge field in this case will be U{3) — > SU{3) x U{1). Let 
us again consider the large 7 solution for the fast rolling action. It is convenient to define 

^ (4.47) 



\/l - 2XR^ 

which allows us to write the energy and pressure in (14.460 as explicit functions of 7. 
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However the exact solution for the speed of sound can be written as a function of 7 



^2 _ ^7^ - 2\ / f(l + 8y)-4(l + 6y) \ 

which, unlike the other solutions studied so far, ahows for the possibihty that the sound speed 
is zero. This is the case if either of the fohowing critical conditions are satisfied: 

7^ = 2 or .2^4(1 + 6Y) 

7c or, Tc (1 + 8Y) • ^ ^ ^ 

The first condition corresponds to (p'^/h^ = T3/4. The second condition is a little more difficult 
to deal with due to the potential (/)-dependence of Y. Now we can consider the two simplifying 
limits, (i) In the limit where y ^ 0, we see that the constraint on the velocity becomes 
02//j4 ^ 3T3/8, while (n) in the converse limit (where Y is dominant), we see that (p'^/h^ = T3/3. 
All of these conditions are allowed because they satisfy the causality constraint on the velocity. 
We know that fluctuation modes exit the horizon at the reduced scale kCg = aH in these 
models, so a zero speed of sound tells us that the modes will never exit the horizon. In order to 
consider an inflationary epoch, we have to ensure the velocity is either much smaller than either 
of the critical bounds (corresponding to non-relativistic motion), or much higher corresponding 
to ultra relativistic motion. Thus unlike the case of = 2, we are lead to selecting a specific 
velocity range. Even for Y ~ 0{1), we find that Cg rapidly tends towards unity as in normal 
models of scalar field infiation. 

In order to consider infiationary solution, we start with the continuity equation. Taking the 
large velocity, we find the general result 



„ -(±1)M^H'J2(1 + 2Y) 

where we have made use of the fact that (p = ±\/2fTsh'^ + .. in this limit. The sign ambiguity here 
can be resolved by demanding 7 to be positive. Since we are interested in solutions where H' < 0, 
we take the + sign in the definition of the velocity. Substituting our expression back into the 
Friedmann equation, using the formula (I4.46|) for the energy density, yields the Hamilton-Jacobi 
equation 

^_ = y(^) _ /,4 (4.51) 

which can again be integrated to solve for H once we specify the background potential. The 
level of non gaussianities arising from this action can be parameterised by 

162 / 10(1 + 8y) n\ 

which clearly has non-trivial dependence on the inflaton fleld for the mass gap backgrounds 
(since Y is constant for the Ads solutions). Let us explore the possible solution space here. For 
the Ads case we find that 

Y = — (4.53) 
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and so can be small with appropriate tuning of the fluxes and the string coupling. If we assume 
y << 1 then we see that the non-gaussianities are (up to 0{Y'^) terms - and dropping the 
constant piece) 

fNL ~ 0.3772 + . . . , (4.54) 
whilst if we assume that Y is large (corresponding to large fluxes) we flnd the following 

fNL ~ 2.997^ + O (^y^ . (4.55) 

The latter condition is much larger than anything encountered before, and severely restricts the 
relativistic approximation we have been making. In fact if we have Y ^ 0(1) we find a similar 
condition. However for small Y we see that the non-gaussianities are roughly the same as in the 
previous sections, and would appear to be the more favourable regime for inflation. This tells 
us that we require Qs « 1/ {MK) . 



4.3.1 Inflation in AdS type backgrounds 



It is generically difficult to find inflationary solution for AdS backgrounds. To proceed with our 
small (j), but large gamma solution, we again turn to a Taylor series approach to the Hubble 
parameter. Let us take the same form for the expansion as in the last section, with a similar 
expression for the inflaton potential. Again we flnd that there is no linear dependence in this 
limit, but the non zero coefflcients can be seen to be 



^0 = J^, (4.56) 



H2 
H3 



M2 ' 

V2T3 
'4M2//o' 



The conditions for inflation are basically the same as in the previous section, which is to be 
expected since we are assuming that inflation is essentially driven by the constant part of the 
Hubble parameter. The slow roll parameter is shifted only slightly by the extra brane because 
the velocity in this case is increased by an extra factor of \/2- Solving for the inflaton at the 
end of inflation, we find 



2 u2 \ 1/3 







(4.57) 



mm 

where we have absorbed the minus sign into the definition of \H2\ to make the solution manifestly 
positive. This is only slightly different from that obtained in the N = 2 case. By integrating the 
Hubble term, we can invert again the resulting expression to obtain the inflaton as a function of 
the number of e-foldings. The result is the same as for N = 2 expect now the tension is doubled. 
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Finally we obtain e as a function of the number of e-foldings 



1+ A.. (1^1)" I (4.58) 

which represents only a slight numerical shift with regard to the expression in the previous 
section for N = 2. We can once again calculate the relevant signals for this model, and the 
analysis proceeds much as in the case of = 2, except that we are now forced to restrict 
ourselves to the 7 >> 1 solution. The tensor amplitude at leading order becomes 

which is a factor of 3/2 larger than the amplitude in the N = 2 case (|4.25p . Whilst the scalar 
amplitude can again be written solely in terms of the tensor amplitude divided by the parameter 
r. 

The tensor spectral index is relatively suppressed as in the N = 2 case, however for the 
scalar index we find 

which is identical in form to the N = 2 solution in (14.280 under the replacement = 
(4Mp/L^l^)^/'^. The same remarks apply here except the physical constraint is slightly tighter 
than before 

> Ml (4.61) 
3.2 xlO"^ - P ^ ^ 

Using the small Y constraint in order to suppress the non-gaussianities we can write this con- 
straint purely in terms of the potential term 

gs X lO^M^ 

V2 » ^^^4 « 1) (4.62) 

Let us consider the two limiting solutions, bearing in mind that we expect that the Y >> 1 case 
will lead to extremely large non-gaussianities. For small Y we see that the sound speed becomes 



.772 _ 4 

Cs{Y « 1) ~ (4.63) 
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and so if we also assume that 7^ >> 4 then we see that this becomes 1/3. This is unlike all the 
other DBI models studies so far. Interestingly if we take the limit where y >> 1 we also see 
that it drops out of the analysis 



CiY » 1) ~ (4.64) 

7 

and in fact we find that ~ 1 as in standard slow roll models of inflation. Following the same 
procedure as in the N = 2 case we can constrain the potential using the scalar amplitude. We 
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find that the range for Vq is 



Vo ~ 0(10-10 _ 10-9) f Y«l (4.65) 



Ms 

which once again indicates the sensitivity of inflation to the string scale. 
4.3.2 Inflation in mass gap backgrounds. 

Let us now restrict our analysis to the mass gap backgrounds. We can again solve the master 
equation assuming that the constant part of the potential dominates the solution. The result is 

Hi<^) = »HMtanh { 4t±^yV^^ , (4.66) 
^ ' Mp y V8Mph2 V u y , y J 

which should be valid for small values of the inflaton field, and we have left the mass gap warp 
factor as an arbitrary constant. We must ensure that this expression is consistent with our 
demand that the 7 factor is large. This requires us firstly to take the minus sign in the velocity 
term, and secondly to examine the behaviour of the function for small values of the inflaton. 
Differentiating this function and then performing a Taylor series expansion of (|4.50|) for small (j) 
yields the constraint (valid up to terms of 0(0^)) 

where we have introduced the simplifying notation Q = Vo/h'^ — 1. Clearly to satisfy this 
condition, we must require that the term C arising from the boundary condition be very small 
in Planck units. Neglecting the (j)^ terms, we can again use a Taylor series expansion, this time 
for C ~ 0. This corresponds to a specific choice of the boundary conditions for the solution. In 
the leading order, we must satisfy the following condition on the parameter Q: 

(4.68) 

However the fact that we require C to be small has an effect on the amount of inflation we can 
have in this fast rolling regime. To see this, let us calculate the fast rolling parameter e, making 
use of our near relativistic approximation. A short calculation shows that 

For small values of C and the inflaton field, we see that the real part of this function is divergent. 
In fact e is a decreasing function of (j) which suggests that inflation will only be possible once 

the field reaches a critical point given by (pc ~ Y^Arccsch(Y'^) — C, after which we enter 
a phase of eternal inflation which will not end within the bounds set by our theory. Of course 
this may no longer be true once higher order terms are included, and we leave this possibility 
for future work. 
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4.3.3 Non-relativistic limits 



In this subsection, we examine the non-relativistic motion of the branes and compare with the 
resuhs from the previous sections. /,From the continuity equation and the definition of the 
energy momentum tensor, we find the inflaton velocity 



2M^,H'y/TTW 

(t>= ^ • 

^ 3(1 + 4y) 

Upon substitution of this back into the Friedmann equation, we obtain 



(4.70) 



2T3 vi + 2y \ ) 2 

where we have introduced another function Z3: 

(3 + 16y)(l + 2y) 



1 + + ,4 + o (^('A) - /i'), (4.71) 



(4.72) 



(l + 4y)2 

The above equation ()4.7ip is difficult to solve analytically for either background. So we resort 
to the usual trick of Taylor expanding the Hubble term for a given potential. 

Let us consider the mass gap backgrounds. The simplest analytic solutions are obtained 
when we keep only terms up to quadratic in the potential and Hubble parameter. It is easy to 
see that the coefficient Hi is imaginary and so we drop it from the analysis. The results for the 
remaining components are 



The general result for the e equation reduces to 

8Mp2g20^ yi + 2Y 
^ 3Hl (1 + 4y) ' 



(4.75) 



for all backgrounds, where y is a function of the inflaton for the mass gap solutions. For small 
values of <j), we can expand this and obtain a value for the field at the end of infiation. As a 
result we can write the slow roll parameter as an explicit function of the number of e-foldings: 

e~ef^% (4.76) 

where we have used the previous definition of /3. For any inflation to occur we must have 
/3 < ^ H2 < therefore we must again choose the minus sign in the expression above. 

Following the same line of reasoning as in the N = 2 case we see that the constraint on the 
potential contributions can be written as follows 
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where we have neglected higher order contributions in /i^. This is smaher than the constraint in 
the two-brane solution due to the additional 'mass' coming from the extra brane. The inflaton is 
effectively weighted by this contribution thus making it roll more slowly. Of course this analysis 
is only valid for small velocities, which is more problematic in this instance as there are zeros 
for the speed of sound function which destroys any hope of obtaining an inflationary solution. 



We can also obtain a simple analytic solution if the constant parts of the potential and the 
Hubble parameter are the dominant contribution. In this case, we find 



Ho = J^(vo-h^ + ^), (4.78) 



, M2 V 3 

where the warp factor contribution is subdominant. This implies that the velocity of the inflaton 
will be zero, as can be noted from the Hamilton-Jacobi expression. 



As for the AdS type backgrounds, we again cannot obtain simple analytic solutions when we 
keep quartic terms in the Priedmann equation. So again we restrict our analysis to the purely 
quadratic pieces. As in the other cases, the linear term in H must vanish for consistency, and 
so the physical solutions are 



-0 - m 



The inflation in this limit is parameterised by the slow roll term 



"^^"H 3(i + 4y) J' 

where we have reintroduced the parameter 6 as in the N = 2 section, which is once again given 

by 

_ MlH2 
Ho ■ 

The validity of the expression is determined by the background fluxes manifest in the Y terms. 
Again we find a similar bound on V2 as in the two-brane case, namely 

V2 > ^^F{Y) (4.82) 
where the function F{Y) now ranges between F(Y) = 1 — > 2\/2Y. 



In both cases we see that in order to satisfy the observed scalar curvature bound, we require 

which can again be used to constrain the maximal value of V2. Once again we see that the 
inflationary scale is sensitive to the magnitude of the string scale. 
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5 Discussion 



In this paper we have examined the evolution of the IR DBI inflation model in the context 
of the non-Abelian world-volume theory. The extra stringy degrees of freedom present in the 
theory serve to induce a so-called 'fuzzy' potential term which makes the form of the action 
significantly different from that of a single brane when we switch to the finite N theory. We 
have seen that the large N limit of the solution is remarkably similar to the single brane case. 
This is not entirely unexpected since this has a dual description in terms of a single classical 
object due to the Myers effect [8]. In fact we expect that our model to be the field theory dual 
of a supergravity solution consisting of a single Z)5-brane wrapping a non-trivial two-cycle with 
units of magnetic flux on the world volume [21]. However in order for this configuration to be 
stable we will generally have to ensure that there is additional worldvolume electric flux on the 
wrapped brane. From our field theory perspective, consistency of the two descriptions requires 
us to turn on a U(l) electric field on the coincident Z)3-branes, thus breaking the symemtry 
group U{N) SU{N) x U{1). It would be useful to investigate this theory and see if it does 
match up with our non-Abelian one at large N, not only because of the intrinsic interest of how 
gauge fields affect cosmological behaviour, but because this may also leads to new understanding 
of Yang- Mills theory via gauge / gravity duality. 

Within the remit of our assumptions, we have tried to argue that inflation is reasonably 
generic. By this we mean that we can generate at least 60 e-folds with relative ease due to large 
N suppression. However for the Ads type solutions we see that the scalar index cannot satisfy 
observational bounds and is thus immediately ruled out. The mass gap case is a little more 
subtle, because it can be a viable model provided we fine tune the solution. 

The most intriguing results concern the analysis of the models where N is finite. Due to the 
highly nonlinear nature of these theories its is not clear what general conclusions we can draw. 
However we studied in detail the two examples of A^ = 2 and A = 3. (It is certainly possible 
to analyse larger values of A^, but the form of the action becomes more and more complicated). 
Already in these examples, we find dramatically different behaviour from the cosmology of a 
single brane. Most notable is that the speed of sound becomes a far more complicated function 
of the inflaton and we are restricted to certain regions of the moduli space of solutions. For the 
A^ = 3 case we noted that the speed of sound actually has zeros in its function. The physical 
explanation for this is not immediately obvious, but could lead to important physical predictions 
which distinguishes this model from all other string inspired inflationary scenarios. An obvious 
question is whether this phenomenon simply arises from the symmetrised trace conjecture, or 
whether it is an inherent feature of DBI models with an odd (not including l)number of branes. 

Another important result concerning this particular solution is that small values of QsKM 
are prefered which can be satisfied with a small string coupling constant and small fluxes. If we 
can trust string perturbation theory then we can find a decent inflationary model. However for 
large fluxes this leads to large non-gaussianities, which may or may not be a good test of the 
model. Only observations from the Planck mission will resolve this issue. Again one may ask 
whether this is an artifact of having an odd number of branes, as the N = 2 case appears to 
have no such restriction. 
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We also found that the scalar spectral index can never be greater than unity in these models. 
In the case of mass gap backgrounds we found that eternal inflation appears to be a common 
feature of both N = 2 and = 3 models. 



Most of our results required us to make the simplifying assumption of 7 >> 1 in the non- 
Abelian DBI theory, in order to be able to carry out analytical computations. This assumption 
was also used in the Abelian DBI cosmology in [3]. However we were also able to consider the 
non-relativistic regime where 7 can be expanded in a power series about unity. In this limit we 
see that inflation does indeed end in both types of background. However the speed of sound is 
essentially unity and therefore doesn't provide us with a nice falsifiable prediction - since we can 
always tune the parameters to allow for enough inflation and suitably low scalar amplitude. Of 
course this will inevitably change if we consider a full string compactification. A more general 
analysis would probe the intermediate velocity spectrum, which we hope will be pursued in the 
future. 



Our model differs significantly in some respects from a single probe brane, despite sharing 
many similarities. The main difference is that now the world- volume theory is playing the role 
of the universe, rather than just the inflaton sector in the standard IR approach to inflation. In 
fact this makes the model a hybrid between that of DBI inflation [3] and Mirage Cosmology [12]. 
As a consequence one may ask about reheating in our model 0. We cannot simply assume that 
there is another throat which contains the standard model since this will not couple to the 
inflaton sector through open string modes. One potential solution to this problem is to assume 
the existence of angled branes in another throat. As our non- Abelian brane configuration travels 
up the throat it will eventually feel a potential being sourced by these branes and will move 
towards them. The standard model can be realised once these branes intersect. Of course this 
implicitly assumes that we have enough branes in our model to allow for a sufficiently large 
gauge group to encompass the standard SU{3) x SU{2) x U{1) gauge symmetry of the standard 
model. An interesting point is that we could actually begin with a model which has a large 
enough gauge group to contain the Standard Model, but still within the realms of the finite N 
expansion. In this case we may find that inflation direclty mixes with the observable sector, 
and a more comprehensive analysis using [23] must be undertaken if we are to place the model 
on a firm experimental footing. This is an interesting, but challenging subject which is left for 
future study in the context of finite N solutions. We remark that due to the nature of our 
approximations, we anticipate that there is universal behaviour for any value of N for some 
of the observables. For example we conjecture that the scalar index imposes the approximate 
WMAP bound 

> (5.1) 



MK}^MfW_ , 
9s (N + l)- P 



which will clearly be harder to satisfy for larger values of N without tuning. Once again note 
the presence of the infiaton 'mass' term above. This expression is an order of magnitude ap- 
proximation, as we have anticipated that the five-dimensional volume is 0{1)tt^. 



As regards future directions, it would also be useful to embed our model into a more realistic 
compactification along the lines of [22,24], by including the position of the D3-branes in the 
Kahler potential and the non-perturbative superpotential. For the case of finite N, one would 

^°Here is a partial list of papers which deal with reheating [11]. 
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hope that the closed string corrections would be small and computable in a certain region of 
the moduli space (see [25] for the open string computation). The relevant scales will now be 
set by the specific choice of compactification. Since the metric will now also explicitly contain 
couplings to the Kahler moduli, it seems likely that DBI inflation will be more complicated in 
a fully compactifled theory due to the fact that both the radion and the modulus will appear 
inside the DBI action. It may still be possible to obtain inflationary solutions in the single brane 
case, although moduli stabilisation effects may reverse the effects or the warping. In any event, 
one would hope that this could be tackled as a future problem. 

It is also possible to extend the work initiated here to include world-volume Abelian or 
non-Abelian gauge flelds. It is an experimental fact that we see weak magnetic fields at large 
scales in our universe today. Although including a non zero B field complicates the inflationary 
analysis, it will certainly play an important role in a reheating phase. It would also be useful to 
consider an extension of [6] to the case of finite N. Finally, another avenue of possibility would 
be to consider the coupling of open string tachyons in the model, which have been shown to lead 
to interesting inflationary solutions [26]. 
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6 Appendix 



In this section we explicitly calculate the relevant perturbation amplitudes for the non-Abelian 
action. The deflnitions of the parameters in this section differ from those in other sections in 
order to simplify the calculations as much as possible - and we use units where Mp = 1. The 
action in the general case can is a non-linear function of the inflaton fleld and its time derivative, 
therefore it can be written in the following form - consistent with the general prescription 
described in [27, 28] 

^ |+p(<^,X) , (6.1) 
where 



S 



P 



1 - 2h~^{(p)T^^X^l + C~^h' 



(6.2) 



with X = <j)^/2 and C = X^CT^/A. We wfll expUcitly consider the case of large N in this 
appendix, however it is straightforward to show that the derived results also apply for the finite 
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N case. The background equations following from this are 



3H^ = 2Xp,x -P = P, 



1 



3(a'^#,x)' -P,^ = 0. 



Note that the energy density p is given here by 



p = NT3 
where the fuzzy potential is 



1 - 2h~^{cl))T^^X 



h\cP) + V{<P) 



w{ci)) = Vi + c-i/i- 



(6.3) 



(6.4) 



(6.5) 



We consider the following general perturbed metric about a FRW background 



-(1 + 2A)At^ + 2adiBdx'dt 

[(1 + 2'ilj)6ij + 2dijE + 2hij] dx'dx^ 



(6.6) 



where di represents the spatial partial derivative d/dx^ and dij = VjVj — (l/3)(5jjV^. Here A, 
B, ijj and E denote scalar metric perturbations, whereas hij represents tensor perturbations. 
Defining the so-called comoving perturbation 



7^ = ^- 



(6.7) 



the Fourier modes of curvature perturbations satisfy the following expression 



where 



z 



2 a^(P^ip,x + 2Xp^xx] 



V = zTZ , 

P,x 

P,x 



4 



P,x 



px + 2Xp, 



XX 



(6.9) 



Note that A; is a comoving wavenumber and a prime represents a derivative with respect to 
a conformal time t = J a~^dt. If the variable z has a time-dependence z oc \t\'^, one has 
z"/z = 'Js/t'^ with 75 = q{q — 1). As long as c| is a positive constant or a slowly varying 
positive function, the solution for ()6.8p is given by 



ci{k)Hi'J{csk\T\) + C2{k)Hi^Xcsk\T 



(6.10) 



where vs = y^75 + 1 /4 = \q — 1/2 1. The coefficients are chosen to be ci = and C2 = 1 to 
recover positive frequency solutions in a Minkowski vacuum in an asymptotic past. 
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Defining the spectrum of curvature perturbation as V-ji = k^\Tl\'^ /2tt'^ , we obtain 

„2 -Si's / IT\2/ 1 \2 / -p/ , X X 2 / , I IN 3-21^3 



/ffV/ 1 y fr{,.s)Y (k\r\ 

\27rJ yaH\T\J \t{3/2) J 
- (6.11) 



Tlie spectral index of the power spectrum is 



- 1 = 3 - 21/5 = 3 - \/475 + 1 , (6.12) 

which means that the scale-invariant spectrum corresponds to i^s = 3/2. About the de-Sitter 
background with |r| = 1/aH, the amplitude of the curvature perturbation is given by 

1 / X 2 



A's • (6.13) 

P,xcs \2^4>J 

The tensor perturbations satisfy the same equation as in the case of standard slow-roll inflation. 
Taking into account polarization states of tensor modes, The power spectrum is given by 



3-2I/T 

/ // \ / I \ / iiryr i \ i t\.\ i \ \ 

Vt 



/ffv/ 1 y (n.r)y (k\r\ 

\27rJ \aH\T\J \t{3/2) J [ 2 
A'r(^X~'^\ (6.14) 



where ut = ^J^t + 1/4 with a"/« = "It/t"^ ■ Hence about the de-Sitter background the ampli- 
tude of the tensor perturbation is 

Al^s[l^)\ (6.15) 

The spectral index of the power spectrum is 

TOT = 3 - 2z^T = 3 - v^47T + 1 . (6.16) 

The tensor to scalar ratio is 

r = = ^jpP,xcs ■ (6.17) 

To study the running of the spectral indices we find it convenient to introduce the following 
parameters: 

H 4> F , , 

where F = p^x + 2Xp^xx and ei is the same as the e which we use in the main part of the 
paper. If = 0, we can derive 

z" m (l + ei+e2 + e3)(2 + e2 + e3) 

— = — , ln = -^^ \2 • 6.19) 

z T'^ (1 — e\Y 
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Under the slow-roll approximation |ej| <^ 1, we find that the spectral index of the curvature 
perturbation is given by 

ns-l = -2(2ei + eg + £3) • (6.20) 

Similarly the spectral index of the tensor perturbation is 

UT = -2ei . (6.21) 

By using the background equations we have ei = (j)^p^x / {'^H'^)- This then shows that the tensor 
to scalar ratio (|6.17|) yields 

r = 16eiC5 = -^csUT ■ (6.22) 

Again this is the same expression as in the single brane case, and is a distinctive feature of DBI 
inflation. The WMAP normalisation we will employ in this paper are the following [1, 18] 



As^ = 10- 



~9 



^ < 0.55 (< 0.24 at 0.95 C.L) 



Hs 0.987_Q Q2y 



fNL < 100 (6.23) 
which may differ slightly from normalisation used elsewhere. 
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